Low energy excitations of the t—J model in one and two dimensions 
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We present an exact diagonalization study of the low energy singlet and triplet states for both 
ID and 2D t—J model. A scan of the parameter ratio J/t shows that for most low energy states 
in both ID and 2D the excitation energy takes the form E{t, J) — a ■ t + b ■ J . In ID this is 
the natural consequence of the factorization of the low energy wave functions, i.e. spin-charge 
separation. Examination of the low energy eigenstates in 2D shows that most of these are collective 
modes, which for larger J correspond to a periodic modulation of the hole density. The modulation 
is well reproduced by treating holes as hard-core Bosons with an attractive interaction. 



Deviations from 'generic' Fermi liquid behavior ob- 
served in the normal state of high-temperature supercon- 
ductors have inspired a search for exotic ground states of 
strongly correlated electron systems. Thereby a 2 dimen- 
sional version of the Tomonaga-Luttinger liquid (TLL), 
which is realized in various ID systems has received 
considerable attention. In a TLL the entire low energy 
excitation spectrum consists of collective excitations, re- 
sulting in a 'decay' of the physical electron into pairs of 
such collective modes and consequently rather unusual 
physics. 

Exact diagonalization allows to systematically scan the 
entire low energy spectrum of small clusters and even to 
analyze the wave functions of low energy states and thus 
provides a unique tool to obtain approximation-free, un- 
biased results for 2D strong correlation systems. While 
the effects of the cluster size pose a serious limitation, 
it will be shown below that already quite small systems 
suffice to clearly demonstrate e.g. the characteristic fea- 
tures of the TLL in ID. This suggests to seek analogies 
between rings and planar clusters, and, as will be shown 
below, 2D systems indeed show some similarity with ID. 
In particular, in both cases we find a very regular scaling 
behavior of the lowest excitation energies as well as evi- 
dence for a large number of low energy collective excita- 
tions. The collective charge excitations can be modelled 
by treating holes as hard-core Bosons. 
We begin with a discussion of ID. Figure ^ shows the 
dispersion of the excitation energy AE (i.e the energy 
relative to the ground state energy) for the lowermost 
singlet and triplet states of a 12-site ring with two holes. 
Results for both longer and shorter rings (down to 6 sites) 
are completely consistent with these. In the S=0 sector 
one can identify two 'branches' of states with very differ- 
ent dependence on J: a first branch (connected by dashed 
lines) resembles the lower edge of a particle-hole contin- 



uum, and has its lowest excitation energy at 2kp=5TT/6. 
The excitation energy roughly scales with J, although 
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FIG. 1. (a) Lowest excitation energies for two holes in a 
12-site ring for all possible momenta and different J. The fig- 
ure shows the lowest or lowest two states for each momentum, 
the lines are guide to the eye. 

(b) EMD for the ground state of 2 holes in the 12-site ring 
(dashed line) compared to n{q) for some excited states with 
given momentum k (full lines). The ratio J/t=0.5. 

the energies do not really scale to zero as (which 
however is likely to be a finite-size effect). A sec- 
ond branch (connected by a full line) comes down at 
Akp—2'K/6, and is practically invariant under a change 
of J, i.e. its energy scale is t. The lowest triplet states 
closely resemble the dispersion and scaling properties of 
the singlet J-branch, having their minimum excitation 
energy a.t2kF- The singlet t-branch and triplet J-branch 
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have been discussed by Bares and Blatter based on the 
Bethe-ansatz solution at the supersymmetric point, J=2t 
[|| ; the (near) identity of low energy singlet and triplet 
spectrum is known for the undoped 11? Heisenberg anti- 
ferromagnet, where a 'singlet replica' of the de Cloiseaux- 
Pearson mode has been discussed by Mueller et al. 
This suggests that in the language of the Ogata-Shiba 
wave function 0| the J-branches correspond to excita- 
tions of the spin part of the wave function, the t-branch to 
an excitation of spinless Fermions or, equivalently, hard- 
core Bosons. 

The unusual nature of the low energy excitations can 
be nicely seen in the electron momentum distribution 
(EMD), defined as n{q) = (c^jCgj). This is shown in 
Figure |l|b for some of the low energy states with 52=0. 
The EMDs of the excited states differ appreciably from 
that of the ground state only near ifc^, which is in 
principle what one would expect for particle-hole exci- 
tations in a Fermi liquid. The only excited states, how- 
ever, where the change of n{q) would be consistent with 
the Fermi-liquid picture are the ones with momentum 
k=2kF- Here the decrease of n{q) at g=— 27r/6 and 
the nearly equal increase at (j'=37r/6 seem to indicate 
the shift of an electron between these two momenta, as 
expected for a particle-hole excitation with a momen- 
tum transfer of kp=5'K/6. For the excited states with 
other momenta no such assignment of an electron trans- 
fer with the proper momentum is possible (for example 
the 4:kp singlet state 'looks like' a particle- hole excita- 
tion 27r/6^— 37r/6, which however would give a wrong 
momentum transfer of —5tt/6 rather than 27r/6). These 
results may be understood by recalling that the low en- 
ergy excited states may be thought of as particle-hole 
excitations of either a spinon or a holon, whereas an elec- 
tron near kp corresponds to a combination of low energy 
spinon-holon pairs. Changing the 'occupation numbers' 
of spinous and holons therefore will in turn change the 
n{q) of the electrons near kp. How n{q) changes in de- 
tail, however, depends on the (unknown) 'expansion' of 
the physical electrons in terms of spinon-holon pairs, and 
thus need not at all be consistent with the particle-hole 
picture. 

While the lowermost states correspond to 'pure' excita- 
tions of either spin or charge degrees of freedom, there 
is no reason why there should be no 'mixed' excitations, 
where spin and charge part of the wave function are ex- 
cited simultaneously. This is demonstrated in Figure H, 
which shows the excitation energies, for a larger part of 
the low energy spectrum for both 2 and 4 holes in the 12- 
site ring. Some of the excitation energies (indicated by 
full lines) have the form AE{t, J)=a ■ t+b ■ J. This could 
be understood if one were to assume that for these states 
both, charge and spin degrees of freedom are excited, the 
former having t as their energy scale, the latter having J. 
When two or more of the linear AE vs. J curves inter- 
sect, they sometimes produce pronounced 'hybridization 



gaps'; we believe this is the case e.g. in the singlet sector 
for two holes and momentum 57r/6. 
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FIG. 2. Lowest excitation energies AE for the 12 site ring 
with 2 holes (top row) and 4 holes (bottom row). Spin and 
momentum quantum numbers refer to both panels of one col- 
umn. 

It should be noted that states in one panel are identical in 
all their possible quantum numbers; nonvanishing matrix 
elements of the Hamiltonian are therefore possible, and 
depending on the strength of these matrix elements levels 
may 'repel' each other appreciably. In addition, there ap- 
pear to be deviations from the 'ideal' behavior for large J; 
we believe that this is due to the increasing tendency to- 
wards hole clustering, which in the thermodynamic limit 
is known to lead to phase separation for J larger than 
Jc«3.5 [Q. In the small clusters there is no true phase 
transition, but rather a very slow and gradual crossover 
to the clustered state, which may lead to a 'bending' of 
the AE vs. J curves. The level crossing at J/tf^A, where 
the 4 hole singlet state with fc=57r/6 comes down below 
the ground state, most probably is also related to hole 
clustering. Summarizing the results in ID, the energies 
of most of the low lying eigenstates to a good approxi- 
mation can be written in the form AE{t, J)=a ■ t+b ■ J. 
This can be understood as a consequence of the factoriza- 
tion of the eigenfunctions into a spin part and a spinless 
Fermion part, as established by Ogata and Shiba in the 
limit J/t 0. 

Motivated by these systematics we have performed a scan 
of the J/t dependence of the lowermost excitation ener- 
gies in 2D. In 2D point group symmetry presents a slight 
complication, in that all allowed momenta for both the 
16 and 18 sites clusters are along high symmetry lines; 
the respective eigenstates then can be classified accord- 
ing to their parity under reflection by the high-symmetry 
line. In the following we have only selected those states, 
whose parity equals that of the ground state, so that 
they could in principle be seen in the dynamical charge 
correlation function (DCF). Then, Figure || shows some 
representative data for the J dependence of the lowest 
singlet and triplet excitation energies. As was the case 
for ID some of the excitation energies (indicated by full 
lines) have the form AE{t, J) = a ■ t + b ■ J . As in ID 
there are dense continua (indicated by dashed lines) with 
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energy scale J coming down at small J so that the lin- 
ear branches 'get submerged'. For some of the states the 
Ai? vs. J curve has a negative curvature; we believe 
that this is a manifestation of the stronger tendency to- 
wards hole clustering in 2D. Extending the range of J 
up to extreme values such as J/t=10, where the energy 
is essentially determined by the number of broken bonds, 
the energies of these states approach that of the ground 
state. We therefore believe that the negative curvature 
is a consequence of increasing hole clustering, which in 
the thermodynamic limit leads to phase separation for 
J greater than Jc~>l.5 Since this critical value is 
smaller in 2D than in ID it seems plausible that in 2D 
the excitation energies are influenced stronger. Another 
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FIG. 3. Lowest excitation energies A_E for different 2D 
cluster sizes A*' and electron numbers in 2D. Each row 
corresponds to one momentum, the graphs are labelled by 
electron concentration Ne/N and spin S. 

noteworthy point is that previous studies |^ indicate that 
most of the states which contribute significantly in the 
DCF seem to be 'pure i-branches'. 

Despite some complications the above data show an un- 
deniable analogy between the low energy spectra of 1-D 
and 2D. To find out more about the nature of the 2D 
eigenstates, we proceed to an analysis of their wave func- 
tions. We start with a discussion of the static hole density 
correlation function, defined as 

=^((l-nO(l-n,+R)), 

i 

where Ui is the operator of electron density at site i. The 



expectation value is taken with the respective eigenstate, 
some representative results are given in Figure ^ In the 
ground state (with momentum (0, 0)) for J/t=l, g{R) is 
uniform in 16 sites and even shows a slight decrease at 
small distances in 18 sites. Contrary to the fairly uniform 
density correlation in the ground states, the low energy 
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FIG. 4. Hole-density correlation function g{R) for the low- 
est S—0 eigenstates with 4 holes/hard core bosons and differ- 
ent total momentum k: (a) t—J model, 16 sites, (b) hard core 
bosons, 16 sites, (c) t—J model, 18 sites, (c) hard core bosons, 
18 sites. The radius of the circle at each R is proportional 
to the value of g{R), the ratio J/t=l for the t—J model and 
V/t—1 for the bosons. 

singlet states show pronounced and systematic modu- 
lations of g{R). In particular, for momentum along 
(1, 0) these states seem to correspond to 'dynamical hole 
columns'. This is reminiscent of the results of Prelovsek 
and Zotos |Q, who studied static four-point correlation 
functions and for large enough J/t found indications of 
column-like hole correlations in the ground state. Our 
scan of the low energy spectrum indicates that the re- 
spective states develop quite continuously down to small 
values of J/t (for example the state at A;=(7r,0) has a 
completely linear AE down to J=0.5); in this physical 
parameter regime, however, the column-like structures 
become nearly unobservable. In previous studies [^,^, 
the correlation of holes in the ground state was found to 
be very similar to that of hard-core Bosons. By adding an 
additional nearest-neighbor attractive interaction V be- 
tween hard-core Bosons we could indeed reproduce the 
g{K) of the excited states of the t—J model surprisingly 
well, as shown in Figure ^ We believe that this is an 
indication that the columns of holes represent a sound- 
like excitation which gains additional stabilization by a 
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residual short range attractive interaction between holes. 
We also note that such column-like structures are not 
obtained when we excite a particle-hole pair from the 
projected ideal Fermi gas. 

The collective nature of the low lying states again be- 
comes very obvious in the EMD. More precisely, we con- 
sider the difference An{q) of the EMD for an excited 
state minus the EMD of the ground state; this quan- 
tity reflects the redistribution of electrons in momentum 
space when going over from the ground state to the ex- 
cited state. Figure ^ shows An(q) for some of the low 
energy states of the 16 site cluster. As was the case in 113, 
the changes of n{q) occur predominantly in the region of 
k space where the EMD of the ground state drops 
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FIG. 5. Difference An{q) for various low energy states 
in the 4x4 cluster with 4 holes. For each q an open (full) 
circle indicate an positive (negative) An{q), the radius of the 
circle is proportional to the magnitude. The ratio J/t=0.5, 
the An(fc) refer to the lowest singlet (top row) and triplet 
(bottom row) state with momentum k. 

sharply; this region is indicated by the dashed line in 
Figure ^. Moreover, it is usually impossible to assign 
transitions of a single electron which could explain the 
momentum of the state. Let us for example consider 
the An{q) for the triplet state with k — {tt/2,0). Here 
it might seem that an electron has been moved from 
(tt, 0) (— 7r/2, ±7r/2). However, this is inconsistent 
with the momentum of the state, {tt/2,0), so that this 
state obviously cannot correspond to a single particle- 
hole excitation. We have performed a similar analysis 
for other momenta and point group symmetries, as well 
as for 6 holes in 16 sites for both t—J and Hubbard model 
[ p^ , and found consistency with the particle-hole picture 
only in few cases. The situation is very different at lower 
electron concentration, where Fermi-liquid like particle- 
hole transitions can be clearly identified in the EMD for 
all low energy states pO| . Then, summarizing the numer- 
ical results for 2D, there is first of all an obvious analogy 
with ID. Many of the low lying excitation energies have 
a parameter dependence of the form AE{t, J)— a ■ t+b ■ J, 
which in ID reflects the factorization of the low lying 



eigenf unctions. As was the case in ID most of the low 
energy eigenstates in 2D seem to be collective modes, 
corresponding to excitations of several different particle- 
hole pairs. For large enough J some of these develop 
into 'dynamical hole lines'. This may partly explain the 
analogy with ID, in that a line of holes would 'cut' the 
spin background into two parts in much the same way 
as a single hole would do for a ID system. The low- 
est of these column-like states can be followed from the 
large J limit down to values as small as J«0.5 — 0.3, 
and have relatively small excitation energies. It should 
be noted that our clusters allow only for density modu- 
lations that are 'harmonics' of the cluster boundary and 
that finite-size effects most probably tend to increase ex- 
citation energies. In an infinite system these collective 
modes therefore may well appear at low energies and in 
the physical regime of parameters. 

As for the importance of these collective modes for the 
low energy physics one could envisage different scenar- 
ios. It might be that they simply represent an additional 
complication 'on top of a simple Fermi liquid-like state, 
i.e. a kind of sound waves. A second possibility would be 
that, similar as in ID, the collective modes comprise the 
entire low energy spectrum, leading to the well-known 
'decay' of the physical electron, and thus some 2D ana- 
logue of the Tomonaga-Luttinger liquid. Finally, if one 
were to assume that extremely low energetic 'hole col- 
umn' excitations persist in some way in the physical pa- 
rameter range and that these columns' zero-point fluctu- 
ations dominate the low energy physics, one would arrive 
at the 'meandering domain wall fluid' jl^. However, for 
any of these scenarios, the fact that the low-energy col- 
lective charge excitations also have density correlations 
similar to hard-core Bosons suggest a non-conventional 
metallic state in this system. 
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